Computing eigenvalues and Hermite interpolation for Dirac systems with eigenparameter in boundary conditions

Eigenvalue problems with eigenparameter appearing in the boundary conditions usually have complicated characteristic determinant where zeros cannot be explicitly computed. In this paper we use the derivative sampling theorem ‘Hermite interpolations’ to compute approximate values of the eigenvalues of Dirac systems with eigenvalue parameter in one or two boundary conditions. We use recently derived estimates for the truncation and amplitude errors to compute error bounds. Using computable error bounds, we obtain eigenvalue enclosures. Examples with tables and illustrative figures are given. Also numerical examples, which are given at the end of the paper, give comparisons with the classical sinc-method in Annaby and Tharwat (BIT Numer. Math. 47:699-713, 2007) and explain that the Hermite interpolations method gives remarkably better results. MSC: 34L16; 94A20; 65L15